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Abstract 

The statistical properties of a two dimensional lattice of elastic lines in a random 
medium are studied using the Bethe ansatz. We present a novel mapping of the 
dilute random line lattice onto the weak coupling limit of a pure Bose gas with delta- 
function interactions. Using this mapping, we calculate the cumulants of the free 
energy in the dilute limit exactly. The relation between density and displacement 
correlation functions in the two models is examined and compared with existing 
results from renormalization group and variational ansatze. 
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1 Introduction 



Disorder, in the form of quenched impurities in a sample, is likely to modify 
various measurements. However, from a theoretical point of view, character- 
izing the behavior of disordered systems has proven difficult, making this one 
of the most challenging and controversial fields of statistical mechanics. The 
directed polymer in a random medium (DPRM) is one of the rare cases of a 
disorder-dominated system whose statistical properties can be determined in 
great detail both analytically and numerically [1], and thus serves as a proto- 
type for other such systems. Furthermore, its behavior is related to a plethora 
of other problems in statistical mechanics, such as the kinetic roughening of 
surfaces [2] , the hydrodynamics of the Burger's equations [3] , and surprisingly 
the biological sequence- alignment problem for genes and proteins [4,5]. 
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The DPRM is directly relevant to vortex lines in superconductors with point 
impurities. Most research on the vortex phases in these systems has focused on 
the structural properties of the flux lattice in a sample containing many vor- 
tices. Recently, the statistics of single vortex lines was studied experimentally 
on a two dimensional flux line lattice oriented parallel to a thin micrometer- 
sized film of 2H-NbSe2 [6] . Magnetic response measurements show interesting 
sample dependent fine structure in the constitutive relation B(H), providing 
a fingerprint of the underlying random pinning landscape. More generally, 
the random potential is expected to modify the measurement of various ther- 
modynamic and structural quantities, making it desirable to characterize the 
probability distribution functions (PDFs) for the outcome. At one extreme, 
microscopic quantities such as two-point correlation functions are quite sensi- 
tive to randomness, inducing in some cases complicated multiscaling behavior 
of their PDFs [7] . On the other hand, the free energy and other macroscopic 
quantities are expected to have simpler self-averaging behavior, while non- 
trivial finite-size effects are still present in their corresponding PDFs. 

To date, the PDF of the free energy is known exactly only for a single line, 
the DPRM [8-11]. However, as in superconductors, in most cases elastic lines 
appear at a finite density, leading to novel effects resulting from a competi- 
tion between randomness and line interactions. Similar line lattice structures 
occur in domain walls in the incommensurate phase of charge-density-waves, 
in monolayers adsorbed on crystal surfaces, and in the form of steps on the 
vicinal surface of a crystal [12]. As pointed out by de Gennes in the con- 
text of polymers [13], and later on by Pokrovsky and Talapov for domain 
walls in incommensurate phases [14] , the statistical mechanics of a pure two- 
dimensional line lattice maps onto one-dimensional free fermions if the repul- 
sive line interactions are replaced by non-crossing constraints. Applying the 
same constraint, Kardar used the replica method to show that the line lattice 
in a random potential maps onto one-dimensional U (n) fermions with an at- 
tractive interaction [8]. Calculating the ground state energy of this quantum 
system for integer n by Bethe ansatz, and analytically continuing to n — > 0, 
the quench-averaged free energy was calculated in this reference. Although the 
PDF was not calculated explicitly, a scaling ansatz for the replica free energy 
was proposed based on dimensional arguments. Following up on this mapping, 
there have been attempts [15,16] to calculate the density correlations of the 
line lattice from the excitation spectrum of the U (n) fermions by Bethe ansatz, 
leading to different results. 

In this paper we examine the analogy between the replicated line lattice and 
U (n) fermions further. A considerable difficulty in the replica approach is the 
need for analytical continuation to small n of the results obtained for integer 
n. Analyticity is usually assumed in replica theories, since it cannot be proved 
in most cases. Here we develop a new method for analytical continuation in 
the line lattice problem which is based on Euler's gamma function, the natural 
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continuation of the factorial to real valued arguments. This technique provides 
an analytic expression for the replica free energy of the line lattice at small n 
in terms of an integral equation, which in general has to be solved numerically. 
Interestingly, in the dilute limit of the line lattice, the analytically continued 
Bethe ansatz equations are identical to that of the ID Bose gas with delta- 
function repulsions. This novel mapping allows us to obtain exact expressions 
for the free energy cumulants of the line lattice, and to reveal new relations 
between density correlations in the ID Bose gas and the replica theory of the 
line lattice, respectively. The scaling predictions for the replica free energy 
proposed in [8], and used in [17] in a wider context of disordered systems, are 
strikingly confirmed. 

The outline of the paper is as follows. Using the replica method and the 
analogy between elastic lines and world-line of quantum particles, we first 
show that the random line system is described by a quantum problem of 
"colored" fermions, which interact upon contact. This analogy is explained in 
detail in Section 2. Exact solutions for the ground state energy of this colored 
fermion model in the thermodynamic limit are then obtained via the series 
of nested Bethe ansatze introduced by Yang [18]. For an integer number n of 
colors, the resulting sets of Bethe ansatz equations can be reduced to a single 
set of equations by considering n-clusters of fermions, as showed by Takahashi 
[19] and Kardar [8]. These integer n equations are derived in Subsection 3.1. 
In the following Subsection 3.2, we develop a new technique to continue the 
discrete Bethe ansatz equations directly to continuous n. We find that for 
small n, and in the limit of low line density, the complex n-color fermion 
problem can be mapped onto the well-studied one-dimensional Bose gas with 
delta-function repulsions [20] (see Section 4). It turns out that the Bose gas 
interaction strength is proportional to the number n of replicas, thus mapping 
the interesting regime n < 1 into the weak coupling limit of the Bose gas. 
Therefore, the problem of calculating disorder averaged cumulants of the free 
energy of the line system becomes equivalent to a perturbative expansion of 
the ground state energy of the one-dimensional Bose gas in the weak coupling 
limit. The exact results for the cumulants of the free energy are presented in 
Section 5. Finally, in Section 6 we apply our novel mapping to the question 
of what can be learned about the density- density correlation function of the 
line lattice from the known results for corresponding correlation functions of 
the ID Bose gas. 



2 Mapping fluctuating lines to fermions 

Configurations of directed elastic lines can be regarded as the world lines of 
quantum particles, where the coordinate along the line plays the role of time. 
To develop a concrete picture of this analogy, we consider the canonical parti- 
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tion sum Zn of a 1 + 1 dimensional lattice of N interacting lines in a random 
environment of lateral length L. To simplify the sum over all configurations of 
the lines, we assume periodic boundary conditions in the direction parallel to 
the lines. This is justified since the statistical mechanics of the line lattice is 
insensitive to the particular choice of boundary conditions in the limit where L 
is much larger then the spacing between line collisions [21]. Now the partition 
sum can be written as the path integral 

Z N N r rxj(L)=r P . 
W.y)] = wEII dr H J Vx 3 {y)^{-H\V{x,y)}/T}, (1) 

iv ! p j=l J Jxj(0)=rj 

where we have summed over all permutations P and have used Feynman's 
definition of the path-integral measure T>Xj(y) [22]. Here Zq denotes the par- 
tition function of a single line in a pure system with only one end of the line 
fixecQ. The Hamiltonian in Eq. (1) is given by 

which is a functional of the random potential V(x, y). The first term measures 
the elastic energy of individual lines with line tension g. The lines interact via 
a short ranged repulsive pair potential U(x), which is local in y. Effectively, 
a line interacts only with the lines in the same constant y cross section. This 
assumption is valid if the line coordinates Xj(y) vary slowly with y. The last 
term couples the potential from randomly distributed impurities to the local 
line positions, and makes the Hamiltonian "time" dependent. Without this last 
term, the expression for Zn is formally identical to the Feynman path integral 
for the canonical partition sum of an interacting ID Bose gas [21]. But there 
is an additional restriction for the configurations of elastic lines compared to 
the boson world lines. Usually the overlap of lines is associated with a high 
energy cost and, therefore, we treat the lines as non-crossing, corresponding 
to an infinite potential U(x) at x = 0. As a consequence, only the identity 
permutation contributes to the sum in Eq. (1). 

The free energy F[V(x,y)} = — TlnZjv is itself a random variable. To ob- 
tain the probability distribution function (PDF) of the free energy, we use 
the replica method. Moments of the partition function are obtained by intro- 
ducing n replicas of the original system, described by line coordinates Xj tCe (y), 
a — 1, ... ,n, and then averaging over the random potential. Assuming that 
V(x, y) is Gaussian distributed with zero mean and spatially uncorrelated with 



3 Zq depends on the lattice constants of the discretization of the path integral. The 
partition function of the interacting random line system will be expressed relative 
to Z [21]. 
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variance 



[V(x, y)V(x', y')] = A5(x - x')5(y - y'), (3) 

where [. . . ] denotes impurity averaging, the average over the random potential 
yields 

™ = T^ivT II II / dr ^ / ' ' V*:AV) ex P {-HJT} , (4) 

\ ly -J a =lj=l J J Xj, a (0)=r jiC , 

with the replicated Hamiltonian 



Hn = [ dy \± £ | fe) + EE ^(v) - *;,«(!/)) 

J0 [ a =lj=l z \ u i/ / a=li<j 

A n N ) 

E S(x ita (y) - x jtf) (y)) > . (5) 

ZJ a,/3=lij=l J 

The sum over permutations has been eliminated from Eq. (4) due to the non- 
crossing condition. The impurity averaged moments of are obtained from 
the time-independent Hamiltonian H n at the cost of introducing a new attrac- 
tive interaction proportional to the variance A of the random potential. This 
attraction can be understood as follows. In the original random system the 
lines prefer to wander through regions of favorable values of the random poten- 
tial V(x,y). Due to the impurity average the system becomes translationally 
invariant and there are no longer energetically favored regions of space-time. 
Instead, in calculating nth moments of Z^, the system of n copies of the lines 
gains due to [exp(— V(x, y))] = exp(A/2) an energy A/2T per crossing of two 
lines (of different replicas), thus leading to an attraction. 



The time-independence of the replicated Hamiltonian, and the choice of peri- 
odic boundary conditions, allows us to write the partition function in terms of 
the symmetric statistical density matrix p sy m( x i,i, • • • , x N,n] x [ t i, • • • , x N,ni L) 
as an integral over all boundary positions of the lines at y — 0, 

N n 

L). (6) 

j=la=l 



The density matrix itself can be written as an imaginary time Feynman path 
integral for the world lines of Nn particles with fixed boundary positions 
at y = and x' ±1 , . . . ,x' Nn at y = L [22]. The density matrix 
can alternatively be expressed in coordinate representation as 
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Psym . 1 ; • • • j %N,ni •£% 1 ' ' ' ' i N ni -^^) 

= E e-^ E ^i(x hl} . . . ,x N ^*{x' . . . ,x' N ), (7) 



with the mapping /3 qm — > (3 C \L = L/T, h — > = T between quantum me- 
chanical and classical parameters^. The sum runs over all symmetric states 
with energy eigenvalues Ei and eigenfunctions ipi of the 1-dimensional quantum- 
Hamiltonian 



r-p2 n N q2 n ^ n N 

H = —zz E E jnr- + E E u fa - x i) - yp E E - ^) 

a=l j=l ux j,a a=li<j Zi a,/3=li,j=l 



From a classical statistical mechanics point of view, the operator exp(—H/T) 
is the transfer matrix corresponding to the partition function [Z^]. 

Using the representation of the density matrix in Eq. (7), we can perform the 
integrations in Eq. (6) easily. Taking the thermodynamic limit /5 qm ~ L — > oo, 
the free energy F n of the replicated system is given by the ground state energy 
Eq of the quantum system in Eq. (8), 

F n = -T ln[Z™ ] = -TnN ln(Z ) + E L. (9) 

Therefore, the statistical properties of the line lattice in the thermodynamic 
limit L — > oo are dominated by the ground state wave function ipo, which 
must be of bosonic symmetry [22] . To obtain the disorder averaged free energy 
[F] of the original line system within the replica method usually the identity 
[InZjv] = lim n _»o([^j\r] — l)/ n * s used. However, calculating the average [Z%\ 
for all values n > by analytic continuation of E Q (n) yields much additional 
information about the PDF of the free energy beyond [F]. Since [Z^] is the 
characteristic function of the random variable [InZ], expansion with respect 
to n shows that 

where [F^\ c is the jth cumulant of the free energy. This identification relies 
on the possibility of an analytic expansion around n — > 0. That this is indeed 
the case for the line system considered here will be further discussed below. 

Determining the PDF of the free energy is thus reduced to calculating the 
ground state energy of the quantum-mechanical system in Eq. (8). Before we 
can proceed with the calculation of Eq{ji), we have to specify the pair potential 
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U(x). We will consider the limit where the characteristic length scale of the 
short ranged pair potential goes to zero, leaving just a non-crossing condition 
for the lines. This replacement of the interaction potential by a pure geometric 
constraint for the line configurations does not change the physical mechanisms 
involved in the non-trivial characteristics of the line system. The crucial com- 
petition between the elasticity of single lines, the tendency of lines to follow 
energy valleys of the random potential, and the restriction of configurations 
by adjacent lines are preserved. Within the quantum-mechanical description 
of the lines, there is a convenient way to implement the non-crossing condi- 
tion: As noted first by Pokrovsky and Talapov [14] for pure systems, the Pauli 
exclusion principle can be used to represent the elastic lines as world lines 
of fermions, which automatically avoids crossings. We thus look for a ground 
state wave function i/j with a suitable antisymmetry under particle permu- 
tations, which describes Nn particles with just an attractive interaction on 
contact between any two particles. It is interesting to note that we change the 
symmetry of the originally bosonic ground state function by hand by applying 
the exclusion principle. However, the ground state energy we are interested 
in is the same both in the fermion and the impenetrable boson description. 
Introducing the new subscript i — (i — l)n + a instead of (i, a) to number 
consecutively all nN particles, the Hamiltonian (8) can now be written as 



with mass m = g and interaction amplitude Co = A/2T. The explicit distinc- 
tion between lines in the same or different replicas has disappeared from the 
Hamiltonian, but is hidden in the required transformation properties of the 
wave function under particle permutations. The appropriate symmetry of the 
ground state wave function is restricted by two conditions: (i) i/jq must be anti- 
symmetric if two of the N particles in the same replica (color) are interchanged 
(Fermi statistics), (ii) i/j must be symmetric if two of the n different replica 
(colored) particles corresponding to the same line are exchanged since replicas 
are equivalent. Due to its antisymmetry, ip tends to zero if two particles of 
the same color approach each other and the attractive delta-function potential 
becomes inactive. This guarantees automatically that the impurity averaging 
results effectively in an attractive interaction on contact between particles of 
different color only. Therefore, the sum over pairs of lines interacting via a 
delta-function in Eq. (8) yields no contribution if a — /3, resulting in the sum 
over all pairs in Eq. (11). To conclude, the original quantum Hamiltonian (8) 
with U(x) replaced by the non-crossing condition for lines belonging to the 
same replica is equivalent to the Hamiltonian (11), if one solves for a ground 
state wave function of appropriate symmetry. 




2c ^25(xi - Xj), 
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3 The Bethe Ansatz 



In one dimension, a number of interacting quantum systems can be solved 
exactly by the Bethe ansatz [23,24]. Lieb and Liniger were the first to apply 
the Bethe ansatz method to a continuum quantum field theory by solving for 
the ground state energy and the excitation spectrum of the one dimensional 
interacting Bose gas with repulsive delta- function pair potentials [20]. While 
their Hamiltonian is exactly given by Eq. (11), the ground state wave func- 
tion that we are looking for has a more complicated structure since it has to 
describe fermions of n different colors. The case n = 2 corresponds to spin-| 
fermions and was solved some time ago by Gaudin [25] for attractive interac- 
tions, and by Yang [18] for repulsive interactions. Choosing n = 3 colors, we 
obtain a model for nuclear matter with an attractive quark-quark interaction, 
which has been discussed in connection with quantum chromodynamics and 
also studied by Bethe ansatz methods [26]. Since we are interested in an ana- 
lytical continuation of the ground state energy for small n< 1, first we find 
a solution for integer n and then continue it analytically. According to the 
symmetry requirements discussed above, the allowed spatial wave functions 
of the U (n) fermions have to transform like the irreducible representation of 
the symmetric group & n N corresponding to the Young diagram [n N ], which 
consists of iV rows of equal length n. Sutherland has constructed a general 
solution of the eigenvalue problem for Hamiltonian (11) with Co < when the 
wave function transforms like an arbitrary irreducible representation of the 
symmetric group [27]. The method involves a nested series of Bethe ansatze 
and leads to coupled integral equations for momentum densities. The same 
method has been applied later to the attractive case by Takahashi [19] and 
Kardar [8] , leading to considerable simplifications in the solution of the Bethe 
ansatz equations. Before we turn to the interesting regime n <C 1, we explain 
briefly the idea of the method of nested Bethe ansatze, in a simplified version 
for the relevant irreducible representation [n N ] with integer values n > 1. 

3. 1 Bethe Ansatz equations for integer n 

The eigenstates of the Hamiltonian (11) are given by Bethe's hypothesis [23]. 
For < xq 1 < . . . < XQ nN < W, with W the transversal system size, we can 
write the wave function as 



Here Q and P are permutations in & h n, and the numbers [Q, P] form a 
(nN)\ x (nN)\ matrix. In the following we denote the columns of this matrix 




(12) 
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by the vector £p = [-,-P]. Next, we define the permutation operators Vi m , 
which act on the £p. Their effect is to interchange the Qi component with 
the Q m component of the vector £p, i.e., Vi m [Q,P] = [Q',P] where the new 
permutation Q' is defined by Q\ = Q m , Q' m = Qi, and Q'j = Qj for j ^ l,m. 
These operators form a (nN)\ dimensional representation of the symmetric 
group & n N- But this representation is reducible, i.e., for a wave function of a 
given symmetry, many coefficients [Q, P] are identical. Therefore, in the fol- 
lowing we consider implicitly only the subspace of the original £p-space which 
corresponds to the irreducible representation [n N ] of & n N- On this subspace 
the operators Vi m act as matrices of the dimension corresponding to the irre- 
ducible representation. 

The requirement of a continuous wave function, which has cusps wherever two 
of its coordinates are equal in order to compensate the delta-function inter- 
action in Eq. (11), imposes relations between the coefficients [Q,P]. In terms 
of two vectors £p, £p/ corresponding to permutations P and P', respectively, 
which differ by a transposition of % and j, i.e, P\ — % — P' m , P rn = j = P(, and 
Pj = Pj for j 7^ l,m, the relations can be written as 



with the operator 



ip. = Y^ P , (13) 



YlT = x Y . + . Xt x Aj . V lm , (14) 
J A, — Aj — ic Ai — Aj — ic 



and the rescaled interaction strength c = 2mc /h 2 = gA/T 3 . If these relations 
are fulfilled, the wave function (12) is an eigenstate of the Hamiltonian (11) for 
any set of momenta Aj. The allowed momenta are then restricted by imposing 
periodic boundary conditions on the wave function. Expressed in terms of 
corresponding to the identity element / of & n N, the conditions force £/ to be 
simultaneously an eigenvector of the nN matrix equations 

e lA ^£ 7 = X j+lj X j+2j ■ ■ ■ X Nj X lj ■ ■ ■ Xj-ij^i, (15) 

with j = 1, . . . ,nN and Xij({Aj}) = V^Y*- . The Aj which solve these equa- 
tions depend, of course, on the chosen irreducible representation. 



So far, we have considered only the spatial part (12) of the full fermionic wave 
function. The full wave function is given by the product of (12) and the color 
wave function x- The product has to be antisymmetric under simultaneous 
interchange of both position and color of two particles. To guaranty this prop- 
erty, the color wave function x has to transform like the conjugate irreducible 
representation [N n ]. One can think of the basis vectors of this representation 
as all the allowed color sequences of length nN which form by linear com- 
bination the function x- Now, the crucial point of Yang's and Sutherland's 
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method for determining the momenta Aj is to solve the matrix equations in 
color space by regarding the color wave function as a spatial wave function 
describing (n — 1)N distinguishable particles on a discrete cyclic chain with 
N identical vacancies. By making a generalized Bethe ansatz introduced by 
Yang, this new problem can be cast into matrix equations identical in form 
to the original one, but corresponding to the lower dimensional representation 
[iV n_1 ] for the color wave function \- This procedure is then continued until 
the problem has been reduced to one for particles of a single color only. To be 
more concrete, the periodic boundary conditions for x rea d 

NX = x j+ij x j+2j ■ ■ ■ X Nj X ij ■ ■ ■ X j-ijX, (16) 

where 

Hi = <? XlW , (17) 

and X^({Aj}) is obtained from Xij({\j}) by the replacement Vij — > —Vij. 
This change of sign for the permutation operator comes from the fact that in 
general two representation matrices of a permutation P corresponding to the 
original and conjugated representation, respectively, differ by the parity of P 
only. To solve Eqs. (16) we use the generalized Bethe hypothesis of Yang for 
X on the discrete cyclic chain. It consists in the ansatz 

X = E[^mA?^Q 1 )---G(Ag { ) n _ i)jv , 2/Q{ri _ 1)JV ), (18) 
p 

where the integers 1 < uq 1 < . . . < VQ (n _ 1)N < nN denote the coordinates 

of the distinguishable particles, and \ f\ . . . , A^^ is a set of unequal com- 
plex quasi-momenta. P and Q are now elements of @( n _i)jv- The coefficients 
[Q,-P] form now a (n — 1)N x (n — 1)N matrix, the columns of which we 
denote by £' p . Again, we define permutation operators Vi m that act on £p so 
that they interchange Q\ and Q m . By considering a suitable subspace of the 
£p-space, these operators are chosen to form the irreducible representation 
[jyn.-i] Q £ ©.^^ in order to assure that the color wave function has the 
appropriate symmetry. Physically, the number x associated with a particular 
set {yi, . . . , y( n -i)jv} an d permutation Q can be interpreted as the amplitude 
for finding a particular arrangement of colors for the nN particles. Indeed, 
the positions of all the particles of one color, say red, determine the i/i of the 
remaining particles, and the color arrangement of these remaining (n — 1)N 
particles can be identified with the permutation Q. To complete the ansatz, 
we have to specify the function G(X,y), 

The requirement that this ansatz for x is a simultaneous eigenvector of the 
nN matrix equations (16) imposes relations between the £' P , as the original 
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eigenvalue problem for the Hamiltonian (11) does for the £ P . These relations 
can again be written as 



& = y'Up, (20) 



with the operator 



A, (2) - A< 2) 



Y L~ ( 2 ) .(2) . . + x (2) .(2/. . ( 21 ) 

and P, P' defined as above Eq. (13). Finally, we have to impose periodic 
boundary conditions on x- Using the relations (20), we obtain for the vector 
£j the (n — l)N matrix equations 



n t, \ A i A j) c f z C ' — v" y" v" y" d (oo\ 



t^HX? -\j) + c/2 
for I — 1, . . . , (n — 1)N, and with 



x(2) _ » (2) _ . 

*S({A< 2 >}) ee %Y% = *' ( { tC X> 3 ({^}). (23) 

A- — Ay + ic 

Substituting Eq. (23) in Eq. (22), we obtain new matrix equations, which par- 
allel in form precisely the original equations (16) for the color wave function, 
but with operators X'^, in which the quasi-momenta X\ have replaced the 
original momenta \. The new equations read 

— 1 ■ ■ ■ X ( ' n _ 1)iV l X[ l ■ ■ ■ X' l _ ll ^' I , (24) 

with the new set of defined via 

j /" t ; ) ^ ^{ 2) -a( 2 ))-c _^ .(a; 2 )-a j )-c/2 

^ ii i(Ap)-A^ + c M^ (2) -a,) + c/2' [ b) 



Thus Eq. (16) is reduced to the lower dimensional Eq. (24), and the eigenvalue 
fij in Eq. (17), which is obtained from the ansatz (18), is 



( "- 1)Ar ^--Aff) + c/2 
* i(Xj - \ { n)) - c/2 



*= ii r:gr c/ :- (^) 



To reduce Eq. (24) further and to determine the eigenvalue fi[, we apply the 
procedure described above again in the representation [A^ n ~ 2 ] , introducing the 
quasi-momenta A^ 3 \ This process terminates after n — 1 iterations and yields 
Sutherland's Bethe ansatz equations for the Nn(n + l)/2 complex momenta 
with a = 1, . . . ,n and i = 1, . . . , (n + 1 — a)N. Here we have denoted 
the original momenta Aj by A^. 
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Comparing Eq. (17) with Eq. (26), we obtain 



e^ w = J] TTTTil Ik ' T ' with (* = !,-•• ,^). (27) 



(n-l)JV _ \(2) 



i(Af>-Af>)-«/2 



Continuing the iteration, we obtain from Eq. (25) with /xj given by Eq. (26), 
with Xj and \$ replaced by corresponding higher order quasi-momenta, the 
set of equations 



mN ^\( n_m + 2 ) _ ^(™- m + 1 )^ _ c j2 (m-l)N ^(n-m+2) _ ^(n-m+2) 



JJ i(Af _m ^ - Aj n_m+1; ) + c/2 "i W""^ - \\ n - m ^>) + c 



x 




_ A (n-m+3)) + c/2 / i = 1, . . . , (m - l)i\A 



g i(A (—+2)_ A (n- m+ 3) ) _ c/2 . ^ m = n ,...,3 

(28) 



/ 



Finally, the iteration ends at the representation [TV], which corresponds to 
X[j = 1 in Eq. (24). Therefore, this equation is trivially fulfilled with p[ = 1, 
corresponding together with Eq. (25) to 



T— h — "TT^-h rr , with (i = l,...,N). 

(29) 



In the following, we will make a simple ansatz [19,8] for the complex momenta 
\^*\ which solve Eqs. (27)-(29) in the limit W — > oo, leaving only one set of N 
independent equations for the N different real parts of the original momenta 
Aj = A^. To motivate this ansatz, we consider for the time being the case of a 
single line (N = 1) coming in n colors. This corresponds to the representation 
[n] for n bosons without spin. Under the attractive interaction they form a 
bound cluster, which is described by complex momenta X a = ic(n + 1 — 2a)/ 2 
(a = 1,... ,n) forming a so-called "n string" [24]. Switching back to the 
general case of N lines, but now for vanishing disorder strength A so that 
c = 0, the problem becomes that of N free fermions without spin (n = 1) 
represented by the Young diagram [1^]. In this case, the momenta are equally 
spaced along the real axis between — k F and kp = np, where p = N/W is the 
density of particles. With this two limiting cases in mind, it is reasonable to 
make the general ansatz for large W 
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(j = l,...,N), (30) 

« + 7 ) + *<~ m) (31) 

, N; a — 1, . . . ,m + l;m — 1, . . . ,n — 2), 

« + ^) + (32) 
,7V; a = 1,... ,n), 

where the kj are rea/ valued momenta and the correction terms Bj n ~ m ^ will be 
proven to vanish in the limit W — > oo below. In constructing this ansatz we 
start with the quasi-momenta \\ n \ which were introduced in the last iteration 
step of the nested Bethe ansatz. In this last step the original problem of nN 
interacting particles was completely reduced to that of iV colorless cluster 
fermions. Therefore, in analogy to the usual fermionic case discussed above, 
we make the ansatz that the quasi-momenta A correspond to the real valued 
cluster momenta kj. But due to the complex structure of the effective cluster- 
cluster interaction, the momenta kj are, of course, no longer homogeneously 
distributed. To get a physical picture of the internal structure of a single cluster 
of fermions, we think of it as an "n string" composed of n bosons which are 
hold together by the attractive interaction. From the above discussion of the 
N — 1 case we know that the "internal" boson momenta of the cluster are 
equally spaced along the imaginary axis. Therefore, the single cluster wave 
function is localized in space within a range of size 1/c = T 3 /gA. The ansatze 
(30)-(32) show that also for a finite cluster density in the limit W — > oo, the 
overlap of clusters can be neglected, particularly with regard to the internal 
structure of the clusters. Thus each cluster's momenta are simply given by the 
momenta of an unperturbed "n string" plus the collective real valued cluster 
momentum kj. This leads to the nN original momenta A^ of Eq. (32). The 
intermediate quasi-momenta of Eq. (31) can be understood as that of smaller 
sub-clusters which appear during the iterative scheme of the nested Bethe 
ansatze. 



Af = kj + iB™ 
A ( "~ m) -k+iri\ - 

A i=(j-l)(m+l)+a -Kj+lC ^1 

(J=l,... 
^il(j-l)n+a = kj + IC (l - 
(J=l,..- 



We now prove that the ansatze (30)-(32) indeed solve Eqs. (27)-(29) in the 
limit W — > oo, leaving a new set of only N equations to determine the fc/s. In 
substituting the ansatz in Eqs. (27)-(29) we keep the correction terms Bj n ~ m ^ 
only in factors which would become zero otherwise. The resulting equations 
for the correction terms have two important properties: (i) From the first 
set of nN equations [Eq. (27)] one observes that the Bj n ~ m ^ have to vanish 
exponentially fast if W — > oo. (ii) The solutions depend formally on the real 
momenta kj, which are still free parameters. In the following we will show that 
these kj have to fulfill necessarily TV independent equations in order to make 
the equations for the B^~ m ^ consistent and thus solvable. 

Multiplying the n equations of Eq. (27) with the ansatz substituted and i = 
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(i — l)n + a with i fixed and a — 1, . . . , n, we obtain 

ra-1 dC 1 ) _ r( 2 ) V n-1 -/J, k \ I _ _ 

JkinW _ (_-\\n-l TT -"i.jg TT TT H*-* + " C /oo\ 

e - l X J 11 n (i) R ( 2 ) 11 11 _ k \ _ nr - y 66 ) 

In the next step we eliminate all correction terms B^, B^l from this equation 
without calculating them explicitly. Thus we substitute the ansatz into the set 
of Eqs. (28). The cluster momenta kj drop out and we get for to — 3, . . . , n; 
i — 1, . . . ,N, and (3 = 2,... , to — 2, the new set of equations 

^(n— m+2) ^(n— m+1) p(n— m+2) r?(ra— m+2) p(n— m+2) p(n— m+3) 



Tj(n—m+2) p(n— m+1) p(n— m+2) p(n— m+2) p(n— m+2) o( n ~ m +3) ' 

(34) 



In the special case /3 = 1 we have 



Tj(n— m+2) rj(n— m+1) ^(n— m+2) m+3) 



p(n— m+2) p(n— m+1) ^(n— m+2) p(n— m+2) : 



(35) 



and similarly for (5 — to — 1, 

}(n— m+2) 

'i,m—l JJ i,m— 1 _ ■ LJ i,m— 1 ~ J - J i,m—2 



^(n— m+2) m+1) p(n— m+2) p(n— m+2) 



r)( n ~ m+2) p(n— m+1) p(n— m+2) p(n— m+3) ' 
D i,m— 1 — i.m i,m—l ~ -°i,m— 2 



(36) 



Now we take for fixed to G {3, ... , n} and i G {1,... , iV} the product of the 
to — 1 Eqs. (34)-(36). Introducing the following ratio of correction terms, 

m-l g(n~m+2) _ g(n-m+l) 
B ^ m ) = II J'lm+2) „(f-m+l) ' ( 3? ) 

we obtain the simple recursion relation 

Bi(m) = -Bi(m-i). (38) 

From the definition of Bi(m) it is easy to realize that the coefficient of the 
right hand side of Eq. (33) is given by (— l) n ~ l Bi(n) = —Bi(2). Finally, to 
obtain £>j(2) we substitute the ansatze of Eqs. (30)-(32) into Eq. (29) keeping 
again the correction terms only in factors which would become zero otherwise. 
For i — 1, . . . , N, we thus obtain 

13,(2) = ^ = -1. (39) 

This result shows that the correction terms dependent coefficient of Eq. (33) 
has to be one in order to have a consistent set of equations for the B^~ m \ 
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This in turn determines the allowed real parts kj of the original momenta 
Xi = A^ 1 . They have to fulfill Eq. (33), which now reads 

^4n t r !! c+a u=i,---,n). (40) 

¥i 

These are the final Bethe Ansatz equations^ which have to be solved to obtain 
the ground state energy of the quantum Hamiltonian H in Eq. (11) in the 
thermodynamic limit, as 

*=£^=&-"'>" + ££* < 4i) 

i=l 3=1 

In the last equation, we have used Eq. (32). The Bethe ansatz Eq. (40) has 
been analyzed in Ref. [8] to obtain the free energy of the line lattice to leading 
order in the density p = N/W. The central result of this reference is an 
integral equation for integer n, which was obtained from Eq. (40) in the limit 
W, N — > oo with p fixed, by using the fact that [W(kj + i — kj)]~ l becomes a 
continuous function in that limit. By taking the limit n — > 0, and subsequently 
the limit p — > 0, the kernel of the integral equation becomes of Hilbert type, 
thus enabling an analytic solution [28] . To obtain the probability distribution 
function of the free energy, we have to determine the ground state energy as a 
function of n for small but finite arguments. Since extracting this information 
by an exact method from the integer n integral equation of Ref. [8] is not 
obvious, here we start with a direct analytical continuation of the Bethe ansatz 
Eq. (40) to non-integer n. 



3.2 Analytic continuation in n 



In what follows, we show how to evaluate the right hand side (rhs) of Eq. (40) 
for non-integer values of n. The product over the replica index a represents 
simply the extension of a factorial to complex numbers. Therefore, we can 
make use of the recursion relation for the complex gamma function, r(z + 1) = 
zT(z), to obtain for real valued k the relation 



The same equations have been obtained for the case n = 3 in the context of quan- 
tum chromodynamics [26] by a more direct approach which does not use Suther- 
land's nested series of Bethe ansatze. Instead it was a priori assumed that the overlap 
of clusters can be neglected in calculating the momenta in the thermodynamic limit. 
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'iZj 1 ik/c + a _ T(l - ik/c)Y{n + ik/c) 

jJi ifc/c - a ~ ^~ ^ r(l + ifc/c)r(n-ifc/c) 

_ r n n r H fc /c)r(w + tfc/c) 

1 J r(iife/c)r(n - iib/c) ' 1 ' 

which is well defined for all n. Using this representation of the finite product, 
and taking the logarithm of both sides of Eq. (40), the proper extension of the 
Bethe ansatz equations to real valued n reads 



nkjW = —i In 



i=i 



T{n + l {k 3 -k l )lc)T{-i{k,-k l )/c) 
Y(n - iikj - k l )/c)Y(i(k j - fc,)/c) 



(43) 



where some care has to be taken to conserve the property that the ground state 
has a total momentum of zero, i.e., J2f=i kj = 0. To express the argument of 
the logarithm in terms of more elementary functions, we rewrite it as the 
infinite product 

Y(n + ik/c)Y(—ik/c) -pj- m 2 + nm + (k j c) 2 + ink / c 
T(n — ik/c)Y{ik/c) ^ m 2 + nm + (k/c) 2 — ink/c 

Here we have used the product representation 1/Y(z) = ze lz Y\m = i{^-+z/m)e~ z / m 
of the gamma function with 7 the Euler constant [29]. Next we write the last 
expression as a product of exponential factors to compensate the logarithm in 
Eq. (43). Using arctan(z) = ^ln((l +iz)/(l — iz)) we obtain 



Y(n + ik/c)Y(-ik/c) _ ~ 
T(n-ik/c)T(ik/c) "i = V XP 



/ nk/c 

2i arctan 



m 2 + nm + (k/c) 2 



(45) 



This result shows that the rhs of Eq. (43) is indeed a real valued expression, 
and the final Bethe ansatz equation for arbitrary real n becomes 

nk,W = jZg n (^-^, (46) 

with the function 

g n (x) = 2 arctan ( UX — -) . (47) 
^ \m z + nm + x 2 J 

This is the central result of this section. The function g n (x) has no poles for 
real valued x, is bounded between ±71 and converges to ±nn for x — > ±00, see 
Fig. 1. Notice that the jump of 2ir at x = is produced solely by the m = 
term in Eq. (47), whereas the remaining terms are continuous. 

This Bethe ansatz equation is valid for all n and for all line densities p. To 
compare it with results already existing in literature we consider two limiting 
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Fig. 1. Both the function g n ( x ) (solid curve) and its limiting form 
2arctan(n/x) = 7rsgn(x) — 2arctan(x/n) (dashed curve) for 1 « 1, n C 1 show a 
jump of 2tt at x = 0. 

cases. For n = 1 we have g n (x) = 7rsgn(a;), where sgn(x) is the sign function. 
The solution of Eq. (46) then simply consists of momenta kj which are equally 
spaced between —Tip and up for W, N — > 00. This is the free fermion situation 
corresponding to lines which wander due to thermal fluctuations only. Notice 
that with only one color of lines (n = 1) present, the disorder induced inter- 
color interaction is, of course, inactive. The same situation is expected to 
appear in the zero disorder limit c — > 0. Indeed, in this limit g n (x ~ 1/c — > 
00) — > nn sgn(a;) and the n's on both sides of Eq. (46) cancel, leading again 
to the free fermion result. The solution of Eq. (46) in the limit n — > gives 
the disorder averaged free energy of lines in a random environment following 
Eqs. (10) and (41). Expanding g n (x) to first order in n, we get from Eq. (46) 

N 

k 3 W = £ 
1=1 

where we have used the representation 7rcoth(7ra;) = 1 / x +2x J2m=i( m ' 2 '+x 2 )~ 1 
of the hyperbolic cotangent. It is easily seen that this equation corresponds to 
the discrete momenta version of the integral equation derived by a different 
technique for treating the n — > limit in Ref. [8]. In this reference the corre- 
sponding integral equation is solved in the dilute limit p/c — > yielding a free 
energy with steric interaction between lines proportional to p 2 (instead of the 
behavior ~ p 3 due to thermal fluctuations). From Eq. (48), the free energy 
can be obtained at all densities [30] . 

Now consider the thermodynamic limit N, W — > 00 of Eq. (46), while keeping 
the density p = W/N fixed. The employed technique is analogous to that of 
Lieb and Liniger [20], using the fact that the function g(kj) = [W(kj +1 — kj)]^ 1 



kj k\ 



+ ii coth 



7r(kj 



(48) 
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becomes continuous in this limit. The limiting behavior gives the density of 
states in the sense that Wg(k)dk = number of allowed /c-values in the interval 
[k,k+dk]. For a given density p, the boundary (Fermi) momentum K = k^ = 
—hi is determined by the condition 



K 

g{k)dk = p, (49) 
while the ground state energy is given by 

Eo = — n {\ - n 2 )pW + / k 2 Q {k)dk. (50) 

24m v Jl 2m J-k v ' v 1 

Calculating the differences of Eq. (46) between all adjacent momenta kj+i and 
kj and expanding with respect to the difference kj+i — kj, yields an integral 
equation for g(k), 

rK (k — k'\ 

nk = J_ K g n ( — — J g(k')dk', (51) 

for all k G [—K,K\. This integral equation is not amenable to a closed form 
solution because of the complicated form of its kernel g n (x). However, the 
kernel is non-singular, and thus a numerical treatment seems feasible, which 
we leave for a future publication [30]. Instead, in the next Section we will 
show that taking the limit of small p/c leads to an interesting simplification 
which paves the way for a description of the replicated lines in terms of a pure, 
weakly interacting Bose gas. 



4 Mapping to a weakly repulsive Bose gas 



In this section, we develop a novel mapping of the dilute disordered line lat- 
tice onto a weakly interacting pure Bose gas. This mapping then provides a 
valuable tool to obtain exact results for the probability distribution of the 
free energy in the important limits of low density p or high disorder strength 
A ~ c. This is the relevant limit for the critical behavior at the transition 
to a line-free system, e.g., for flux lines expelled at H cl from superconduc- 
tors or for domain walls at the commensurate-incommensurate transitions in 
absorbed layers or charge density waves. 

From Eq. (49) one observes that the cutoff momentum K goes to zero in the 
limit of p — > 0. Therefore the behavior of g n (x) for i < 1 determines the 
solution of the integral Eq. (51) in the limit p/c 1. This limit has to be 
analyzed very carefully since we have to keep the n dependence of g n (x) to 
all orders to obtain a complete expansion of the ground state energy E in n, 
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which is exact to first order in p/c <C 1. Expressing the arctan of the m = 
part of Eq. (47) in terms of its inverse argument, and expanding the remaining 
terms in x, we obtain 

/x\ 00 2n 

9n(x) — 7T sgn(x) — 2 arctan — ) + x — + 0(x 3 ). (52) 

\nj ^ 1 m{m + n) 

Note that we are interested in an expansion of the ground state energy E 
around n = 0. Independent on how small the ratio p/c and therefore the 
argument of g n {x) is, the final expansion in n for fixed p/c tests arbitrarily 
small values of n and thus can render the arctan(x/n) in Eq. (52) of order one. 
Therefore, an expansion of this term is not justified. In contrast, the remaining 
terms with m > 1 in Eq. (52) tend always to zero for both i C 1 and n <C 1 
and can be safely neglected. Thus, retaining the first two terms of Eq. (52) 
only [see Fig. 1] and taking the derivative of Eq. (51) with respect to k, we 
obtain the integral equation 



k 2nc 

k (nc) 2 + (k — k') z 



which yields the complete n dependence of g(k) to first order in p/c exactly. 
The crucial observation is that this equation is identical to that obtained by 
Lieb and Liniger for their exact Bethe ansatz solution of the one-dimensional 
Bose gas with repulsive delta-function interactions [20]. 

To obtain a physical understanding of the relation between the replicated line 
lattice and the Bose gas, we change variables as follows: Define the function 
Qb(k) = g(k/n)/n, the cutoff momentum Kb = nK, the mass mj = nm and 
the rescaled interaction strength q, = n 2 c. With this new variables Eqs. (53), 
(49) read 

2-JT ^ 2 Q b (k')dk' = 2ng b (k)-l, / g b (k)dk = p. (54) 

-Ki C b + yK K ) J— K b 

The bosonic ground state energy, i.e., the second term of Eq. (50) becomes 

h 2 W r K b o 

E , b = / k 2 g b (k)dk. (55) 

2rrib J-K b 

These equations represent the exact ground state (with energy E 0jb ) of a Bose 
gas with repulsive interactions, described by the Hamiltonian 

h 2 N d 2 

Hb = -^T,f02+ 2c °,bY, 6 ( x i- x j)> ( 56 ) 

4111b j = i iK - 

for — > oo with fixed density p of bosons [20]. Due to the definition of 
c = 2mc /ft 2 the amplitude of the interaction is given by c 0) 6 = nc . The cor- 
respondence can be summarized as follows: The exact ground state solution of 
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the quantum Hamiltonian in Eq. (11) for attracting particles, each appearing 
in n colors, is in the dilute limit p/c <C 1, and to all orders in n around n — 0, 
the same as the ground state solution of the Bose gas Hamiltonian of Eq. (56) 
with repulsive interactions. For comparison with the usual mapping between 
elastic lines and world lines of bosons for pure systems, both the correspon- 
dence of parameters for the pure case and the novel analogy of quantities in 
the random system are summarized in Tables 1 and 2. 

Table 1 

The usual correspondence of the parameters of a d — 1-dimensional Bose gas and a 
ci-dimensional lattice of lines without a random potential. The interaction potential 
between vortex lines maps to the Boson pair potential. In d = 2 dimensions, non- 
crossing lines are modeled by hard core bosons or, equivalently, free fermions. 



model 


mass 


Planck's 


inverse 






constant 


temperature 


Bose gas 


m b 


n 




non-random lines 


9 


T 


L 



Table 2 

The correspondence of parameters in the novel analogy between the one-dimensional 
Bose gas and the two-dimensional replica theory of lines. The density of bosons and 
(non-replicated) lines, respectively, is in both descriptions given by p. 



model 


mass 


Planck's 


inverse 


contact 


effective 






constant 


temperature 


interaction 


interaction 


Bose gas 


m b 


n 




CQ,b 

(repulsive) 


Cb 

= 2m h co tb /h 2 


replicated 


nm 


T 


L 


ticq = nA/2T 


2 

n c 


lines 


= ng 






(attractive) 


= n 2 <?A/T 3 



The correspondence between the Bose gas and the replica theory of lines can 
be interpreted physically as follows: As we have seen in Section 3.1, the Bethe 
ansatz solution can be interpreted as n lines of different colors forming a clus- 
ter due to the pair attraction of amplitude Cq. This suggests considering these 
clusters as new particles of mass nm = ng, then appearing as the compo- 
nents of the Bose gas described by Eq. (56). To obtain insight in the effective 
strength of the contact interaction between these composite bosons, consider 
two clusters, each composed of n lines of different colors. Since particles of the 
same color avoid each other like free fermions or, equivalently, bosons with 
hard core repulsion, some inter-cluster interactions are "screened," see Fig. 2. 
There are n\ orderings of particles within a cluster, each ordering providing 
the same number of pair interactions given by 1 + . . . + (n — 1) = \n{n — 1) 
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n=2 — n=3 



om mo 



Fig. 2. The origin of the effective cluster-cluster interaction: For n = 2 there is 
for both configurations only one inter-cluster interaction possible since the black 
particle of the first cluster is blocked by the black particle of the second cluster from 
interacting with the dark gray particle. For n = 3 only two typical configurations 
are shown. It is easy to convince oneself that each of the 3! configurations allows 
for three interactions suggesting 1 + . . . + (n — 1) = \n(n — 1) pair interactions for 
n-particle clusters. 



due to screening. Finally, we have to take into account that in the quantum 
Hamiltonian (11), which corresponds to the replica theory of lines, the inter- 
action sum runs implicitly for two given clusters over all different colors a ^ f3 
instead of running over color pairs a < f3 only. Therefore, the total amplitude 
of the contact interaction between two clusters is given by n(n — 1)cq, which 
yields a repulsion of amplitude — hcq in the limit n <C 1, in agreement with our 
findings (see Table 2). That this naive estimate produces the correct result to 
first order in p/c is presumably related to the fact that the internal structure 
of the clusters does not matter in the dilute limit since the distance between 
clusters is much larger than cluster size ~ 1/c. 

Having established the analogy between the replica theory of lines and the 
Bose gas, one can take advantage of the many results for the one- dimensional 
Bose gas accumulated during the past decades [20,31-33]. In the following, we 
focus on the ground state energy Eq^ as given by Eq. (55) to obtain the free 
energy F n of the replicated line system using Eq. (9). The only dimensionless 
intensive variable in the Bose gas problem is 7 = c b / p = n 2 gA/pT 3 . Thus the 
ground state energy can be written as 



E °' b = ^^ p3e(7) ' (57) 



where e(7) is a dimensionless monotonically increasing function of 7 [20]. Since 
we are interested in the limit n < 1, we have to consider the weak coupling 
limit of the Bose gas. In the limit c b — » 0, the integral equation (54) has the 
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diverging^ solution [20] 



Q b (k) = 2^/*6- fc2 - (58) 

This zero-order result corresponds to the n = limit of g(k), which is a well- 
defined expression according to the mapping between the two models, and 
given by 

lim g(k) = ±Il.y/K*-k>, (59) 
n->Q 2ty gA 

with K = 2^JgAp/T 3 . This result coincides with that obtained in Ref. [8] 
by solving an integral equation which itself is valid for n = only. Using 
Eqs. (54), (55) one gets Kb = 2p r y 1 ^ 2 , e{^) = 7 and, therefore, the non-trivial 
(disorder dependent) part of the disorder averaged free energy of the lines is 
given by 

[F] d = lim ^±L = A (LiL p + A WLt (60 ) 
L i d n^o n 2T \12 T 3 / ' y ' 

using Eqs. (10), (50), (57), and the correspondences of parameters in Table 2. 
Note the difference from the case without random potential, in which the 
non-linear term in [F] is proportional to p 3 . 



5 Distribution of the free energy in the dilute limit 



The preceding Section demonstrates how the lowest order term in the energy 
of the Bose gas yields the first moment or cumulant of the free energy of the 
random line system in the limit p — > 0. As discussed in connection with a 
scaling hypothesis for PDF's of disordered system [34,17], cumulants of the 
free energy are also singular at critical points. To calculate these cumulants 
the complete n dependence of the replica free energy F n is needed. Due to 
the mapping between replicated lines and bosons, F n can be obtained order 
by order around n = by perturbation theory for the Bose gas in the weak 
coupling limit c& = n 2 gA/T 3 <C 1. As realized by Lieb and Liniger [20], and 
as the zero-order solution [see Eq. (58)] suggests, the complete solution of the 
integral equation (54) is highly singular at q, = 0. The physical meaning of 
this singularity in the context of the line system is related to the relevance of 
the random potential. Any randomness, however weak, leads to a glassy state 
that is fundamentally different from the pure case. 

6 The divergence is related to the fact that as q, — ► the kernel of the integral 
equation is a representation for 27r5(k) leading to 2 / KQ\ ) {k) — 1 — > 2ngi ) (k). Therefore, 
Qb(k) has to grow unbounded in this limit. 
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The singular behavior of the integral equation seems to preclude any direct 
analytical way of extracting higher order corrections to the zero-order result. 
Fortunately, for small 7, Bogoliubov's perturbation theory [35] can be used 
to calculate the ground state energy of the Bose gas. One reason for trusting 
Bogoliubov's theory for small 7 is that in the case of a delta-function po- 
tential the ground state energy is given correctly by this theory in the limit 
of high boson density p [20]. Note that the important small n limit maps 
onto the high density limit of the Bose gas as far as the coupling strength 
7 ~ n 2 /p is concerned, although we are studying the low density limit of 
the line system. Indeed, agreement between the exact ground state energy ob- 
tained numerically from the integral equation (54) and Bogoliubov's first-order 
perturbation theory has been confirmed for small 7 by Lieb and Liniger [20]. 
Later Takahashi [31] used the correlated basis function method to calculate 
the second-order correction to the ground state energy, which again coincides 
with the second-order energy in Bogoliubov's perturbation theory [36], and is 
in perfect agreement with a numerical high-accuracy solution of the integral 
equation (54) down to 7 ~ 10~ 5 [37]. Using the analytical results for first and 
second order terms, the next two orders can be obtained by this numerical 
solution, and the ground state energy summarized, in terms of the rescaled 
function introduced in Eq. (57), as 

e(7) = 7 - ^7 3/2 + Q - ±) 7 2 + &57 5/2 + + Otfl\ (61) 
with numerical coefficients [37] 

b 5 = -0.001588, b 6 = -0.000171. (62) 

It is interesting to note that the ground state energy is non-analytic in the 
coupling strength and thus has to be considered as an asymptotic expansion 
for small 7. 

The replica free energy F n of the line system can now be evaluated by applying 
the mapping between bosons and lines. First, we note that Eq. (61) provides 
an explicit proof that F n is indeed an analytic function of n since 7 ~ n 2 . 
While this property of F n is usually assumed in replica theories, it cannot 
be proved in most cases, and had not been checked before for random line 
lattices. In terms of the cluster size 

1 T 3 

'-e = jA' < 63 > 

which sets the crossover length scale from pure to disorder dominated behavior 
for a single line [8], the disorder dependent part of the replica free energy can 
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be written as 



Fn, d = - n 2 )l dP + n{l dP f - ± n 2 (M 3/2 

Z * 2T l 12 37r (64) 



+ 



Q " ^) ^<*P + &5^(/ d p) 1/2 + hn 5 + 0(n 6 ) 1, 



in the dilute limit p <C l d 1 . This exact result for the replica free energy has 
the scaling form 

F n = LWn^- G (n{l d p)- 1 ' 2 ) , (65) 

where Q{x) is a polynomial. This result is in agreement with dimensional 
arguments given in Ref. [8], which motivated a scaling theory for cumulants 
of thermodynamic quantities in random systems [17]. Thus this scaling theory 
is strikingly confirmed by our exact results. 

The cumulants of the free energy of the line lattice at fixed density p can be 
read off from F n . They can be written in the general form 

[F p \c = a p ^AT?- 2 (l d p)^ 2 (66) 



I 2 



with coefficients 



oi = -, cr 2 = — , 0-3 = 7- 4 a 4 = 0.01906, o 5 = -0.01026. (67) 



2 3vr 4 ti 

Consistent with the central limit theorem all cumulants of the free energy 
density vanish in the thermodynamic limit as [f p ] c ~ {LW) 1 ^ P . However, 
for a large but finite system the cumulants are finite and show a non-trivial 
dependence on the density, or on the chemical potential p within a grand 
canonical description. Close to the critical point we have p ~ p — p c and 
[f p ]c ~ (p ~ Pc)^~ p ^ 2 - Therefore, as p — ► p c + the first four cumulants vanish 
continuously at the critical point, the fifth cumulant approaches a constant 
and higher order cumulants show increasing divergence^. 

Of course, since we are discussing finite systems, the singular behavior is cutoff 
when the correlation length approaches the system size. Due to the anisotropy 
of the line system, there are two correlation lengths. The first one is given by 
the mean distance between lines, £j_ = 1/p, whereas the second one is set by 
the mean longitudinal distance between collisions of lines, £y = (gT/A) 1 ^ 2 ^ 2 . 



7 A numerical solution of the integral equation (54) indicates that the coefficient 
by in Eq. (61) is negative and about one order of magnitude smaller than [37]. 
This suggests that the sixth and presumably all higher cumulants have indeed non- 
vanishing contributions ~ WL. 
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The critical regime corresponds to £j_ ^> W, £11 ^> L. We consider a fixed ratio 
of longitudinal and transversal system sizes so that they are in agreement we 
the anisotropic scaling of the system, i.e., L = S^\\(W/^±) 3 ^ 2 , allowing for an 
anisotropy parameter 5. In the critical regime, the moments of the free energy 
density are then given by [f p ] c = a p (A/TW 2 ) p yielding universal relative 
cumulants 

= 2%^, (68) 

which depend only on the geometry of the system via the anisotropy parameter 
5. This result is a consequence of the relevance of disorder, and demonstrates 
the complete destruction of self- averaging at the critical point, since otherwise 
the relative cumulants vanish with increasing system size [34] . Numerical sim- 
ulations for systems of appropriate size ratios should be suitable to test this 
universality. Another interesting property of the PDF of the free energy is its 
asymmetry. Since the third and fifth cumulant are both negative, there is a 
higher probability for the system to be in a state with a free energy which is 
smaller than the average value. This asymmetry gets more pronounced as the 
system gets closer to the critical point. 

Some comments on the above results for the PDF are in order. These com- 
ments are of more general nature, and are given here in order to clarify the 
status of the replica method itself. The replica approach provides a tool to 
calculate integer moments [Z n ] of the partition function. However, Z itself is 
unphysical, and the desired PDF of In Z has to be deduced using analytical 
continuation of the moments [Z n ] to non-integer n, and applying the cumulant 
series representation in Eq. (10). There are two potential problems related to 
this kind of procedure. First, an unique PDF can be deduced from the in- 
teger moments of the PDF only if the [Z n ] grow slower than n\ for large n 
[38]. One might conclude that this condition is not fulfilled for the random 
line system because of [Z n ] ~ exp(— E (n)L/T) and E (n) has a contribution 
proportional to n(l — n 2 ), see Eq. (50). The origin of this contribution can be 
traced back to the energy of individual clusters in the Bethe ansatz. Therefore, 
the same situation appears for the single line system. In this context, it has 
been noted that the result for E (n) given above is not expected to be valid 
for arbitrarily large n [39]. One is lead to this conclusion by the observation 
that the cluster size goes to zero if n tends to infinity, reflecting the tendency 
of attractive bosons to collapse. The issue of an attracting ID Bose gas has 
been discussed briefly by Lieb and Liniger [20] in connection with their Bethe 
ansatz solution. They conclude that for the attractive Bose gas it is not clear 
that the Bethe ansatz wavefunction which includes the n-string solution ex- 
hausts all solutions of the Schrodinger equation. When the cluster size Id goes 
to zero, it seems to be no longer justified to construct the wavefunction as 
a product of only two-paricle functions exp(— \xi — X2\/ld), as done for the 
n-string, since they are extremely localized. The second potential difficulty of 
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the replica method has its origin in the order in which the thermodynamic 
limit and the n — ■> limit are taken. To apply the cumulant series in Eq. (10), 
we must expand the analytically continued function ln[Z n ] around n = for a 
fixed system size. However, in calculating [Z n ] in terms of the quantum system, 
first we consider at fixed n the thermodynamic limit to eliminate all exicted 
states. For a single line in random media the two limits do not commute due 
to the non-trivial sample-to-sample variation of the free energy [39], and the 
free energy cumulants have to be obtained indirectly by a scaling ansatz for 
the replica free energy [17]. For the line lattice considered here, the situation 
is different. Now, the extensive replica free energy F n , obtained by first taking 
the thermodynamic limit and then the limit of small n, has contributions at 
any order in n, and thus the two limits are more likely to commute. 

Finally, we would like to mention that our mapping between bosons and elas- 
tic lines provides also new information about the weak coupling limit of the 
ground state energy of the one-dimensional Bose gas. Since cumulants [F p ] c of 
even order p have to be positive, all coefficients of non-integer powers of 7 in 
the expansion of the ground state energy in Eq. (61) have to be negative, see 
Eq. (10). This is an important constraint on the various approximative meth- 
ods in use for the Bose gas. Recently, different self-consistent approximations 
have been used to calculate the ground state energy based on a so-called local- 
field correction approach [40,41]. In view of the fact that the two approaches 
give different signs for the 7 7 / 2 term in Eq. (61), our sign constraint rules out 
the solution of Ref. [40] from being correct. 



6 Correlation functions 

The preceding study of the free energy of the line lattice shows that disorder 
induces critical behavior which is distinct from the pure case. Moreover, dis- 
order is expected to lead to correlations in the system which are different from 
pure 2D line lattices. In the context of flux lines in high-T c superconductors, 
Fisher predicted a "vortex glass" phase in which disorder locks the flux lines 
into one of many possible metastable states [42] . He showed that the random 
line lattice can be mapped onto the two-dimensional random-field XY (RFXY) 
model with vortices excluded. The correlation functions of the low temperature 
phase of the RFXY model, and the vortex glass phase, and their collective ex- 
citations have been studied extensively during the last decade. For the RFXY 
model, analytical approaches include replica-symmetric renormalization group 
(RG) techniques [43-45], and one-step replica-symmetry breaking variational 
ansatze [46] and RG methods [47,48]. The analytical studies of the random 
line lattice are mainly based on RFXY model results. They consist in scaling 
arguments [49], linear continuum elasticity models [50], RG techniques [51- 
53], and variational methods including replica-symmetry breaking [53] and 
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without use of replicas [54] . However, the results obtained within the different 
approaches are not consistent. Whereas all variational methods predict for 
the continuum displacement field u(x,y) of the lines, correlations of the form 
[((u(x,y) — -u(0,0)) 2 )] ~ ln(a; 2 + y 2 ), RG approaches yield additional correc- 
tions growing as ln 2 (x 2 + y 2 ), which vanish as the non-glassy high-temperature 
phase is approached. A critical comparison of the various results is given in 
Ref. [55]. One should also keep in mind the range of validity of the differ- 
ent approaches. The variational methods are expected to hold throughout the 
whole glassy phase, while the RG results are applicable only close to the glass 
transition point. 

Numerical simulations of glassy line lattices in two dimensions are also in- 
conclusive. For both the variational method predictions [56,57], and the RG 
results [58-61], qualitative agreement was found in simulations. The numerical 
calculations are either hampered by slow glassy dynamics, and/or by a large 
length scale for the crossover to a regime where differences between the two 
analytical methods become significant. More recently, Zeng et al. [62] found 
good quantitative agreement, albeit in the vicinity of the glass transition, be- 
tween the RG results and their numerical studies of the correlation functions 
via a mapping to a discrete dimer model with quenched disorder. 

The analogy between the line lattice and the quantum model of U(n) fermions 
has been used as an alternative to RG for obtaining information about cor- 
relations in the glassy phase [15,16]. References [15] and [16] both use the 
replica-symmetric Bethe ansatz, but employ different techniques in handling 
the n — > limit, and make different assumptions. As a consequence, the two 
approaches yield different results for the asymptotic decay of the line lattice 
density-density correlation function, i.e., [{p(x, y)p(0, 0))] — p 2 ~ (x 2 + y 2 )^ 1 
[15] compared to [(p(x,y)p(0,0))] - p 2 ~ (x 2 + y 2 )~ 1/2 [16]. Although the de- 
cay exponent of the first result is in agreement with the variational ansatze, 
the amplitude calculated in Ref. [15] does not agree with them. Neither result 
shows the logarithmic corrections predicted by the RG approaches. 

In this section we reconsider the displacement and density correlations in the 
glassy phase by exploring the mapping of the random line lattice to the pure 
Bose gas. While the free energy of the former is related to the ground state 
energy of the latter, correlations in the statistical system are related to the 
low energy excitations of the quantum problem. Therefore, we calculate the 
low-energy excitations of the U(n) fermion model in the limit n — > 0. In the 
following, we consider only replica-symmetric or color-less excitations, which 
do not break apart the n-particle clusters corresponding to the bosons. This 
restriction is motivated by the fact that for n > 1 only singlet states are ex- 
pected to be gap-less [63], due to the finite energy associated with breaking 
up a cluster. With this assumption, the excitations above the ground state 
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consist of "cluster-particles" and "cluster- holes." Under the mapping^] of the 
replicated lines or U(n) fermions to the Bose gas, these excitations become 
the usual particle and hole states studied by Lieb [20] . He calculated the dou- 
ble spectrum associated with particle and hole elementary excitations exactly 
by Bethe ansatz. Particle states can be excited with any momentum p, but 
the hole state spectrum exists only for momenta p satisfying the condition 
\p\ < irp. Both spectra have in common a linear dispersion relation at small 
momenta with the same velocity of sound, 

e P M=v s \p\ + 0(p 2 ), (69) 

where v s in the weak coupling limit is given by 




v s = ^ = Jfi**». (70) 
m b V m b 

Since under the mapping both cq^ ~ n and m;, ~ n, the excitation spectra of 
the replicated lines or U(n) fermions remain linear in the n — > limit, and 
the velocity of sound is given by 

^ (71) 

A linear dispersion relation in the n — > limit was also obtained in Refs. [15,16] 
by solving directly a linear integral equation which determines the shifts of the 
Bethe ansatz momenta kj under particle/hole excitations. However, in these 
references different assumptions and approximations were made to go from 
the excitation spectra to the density- density correlation function of the line 
lattice, leading to contradictory results. 

Having the mapping to the pure Bose gas at hand, it is tempting to reana- 
lyze the density-density correlations of the line lattice from a new perspective. 
Below we will show that the mapping gives direct results only for certain cor- 
relations in the replicated system, from which information about correlations 
in the original random system has to be deduced. Therefore, the remaining 
part of the section is guided more by the question of how much we can learn 
about the correlations in the line system from the knowledge of the density 
correlation function of the Bose gas, rather then giving final exact results for 
the line lattice correlations. The density- density correlation function of the ID 
Bose gas at zero temperature has been studied extensively by several methods 
in the past, see, e.g., Ref. [33]. Since the dispersion relation is linear for small 
p, conformal field theory (CFT) can be employed to obtain the asymptotics 



8 That the mapping holds also for the excitation spectrum follows from the fact that 
the response of the Fermi sea (shifts of the discrete wave vectors kj) to excitations 
is completely determined by the phase shifts in an infinite system associated with 
transposing two wave vectors of the Bethe wave function, see, e.g., Ref. [20]. 
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of the correlation function. The decay exponent rjb is given by the conformal 
dimensions, which are determined by the spectra of low-lying excitations. For 
the class of Luttinger liquids, including the Bose gas, the CFT approach is 
equivalent to Haldane's effective harmonic fluid description of ID quantum 
fluids [32]. This description is based on the fact that the elementary exci- 
tations may be regarded as bosons which represent long-wavelength density 
fluctuations. The Hamiltonian for these bosons can be expressed as a sum 
of harmonic oscillators. Therefore, correlation functions can be calculated by 
simple Gaussian averages. In particular, from the Luttinger liquid approach 
and CFT, the asymptotic time-dependent correlation function of the local 
density b(x, t) of the ID Bose gas is known at any coupling strength q, as 
[32,64] 

/, , \w„ ^ 9 Vb x 2 — v 2 t 2 . cos(2irpx) , , . . 

<«*, wo, o» = p> - ^_j-r^ + V + ; >fA^ + hh " (72) 

where the coefficient A is not known exactly and higher harmonics (h.h.) are 
not shown explicitly. The coupling strength enters the exponent r]b only via 
the velocity of sound, which is known perturbatively for small q,, and numeri- 
cally for all q>, from the Bethe ansatz [20]. This exponent also determines the 
amplitude of the homogeneous part of the correlations, and is given by 

Vb = • 73 

m b v s 



We now examine to what extent the knowledge of the density- density corre- 
lation functions of the Bose gas provides information about the random line 
lattice correlations. First consider the correlation function of the replicated 
U(n) fermion system of mean density np. Its local particle density r(x,y) is 
given by the sum over all densities p a (x,y) of particles of a specific color a, 
i.e., 

n 

r(x,y) = ^2Pa(x,y). (74) 

a=l 

Since all replicas or colors are equivalent, we obtain in terms of the density 
fluctuations Sp a (x,y) = p a (x,y) — p, the relation 

(r(x, y)r(0, 0)) = n 2 p 2 + n(5p a {x, y)5p a {0, 0)) + n(n - 1) (8p a (x, y)8pp(0, 0)). 

(75) 

On length scales larger then the cluster size Id, see Eq. (63), the simple rela- 
tion r(x,y) = nb(x,y) between boson and U(n) fermion densities holds, since 
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bosons correspond to n-clusters of fermions. Thus, we obtain 



(b(x, y)b(0, 0)> =P 2 + 1 ((Spa(x, y)S Pa (0, 0)> - (5p a (x, y)5 P/3 (0, 0)» 
+ (5p a (x,y)6pp(0,0)), 

where we have assumed that all off-diagonal correlations (a ^ (5) are identical. 
The disorder and thermally averaged density- density correlation function of 
the original line system can be calculated by taking the n — > limit of one of 
the n equivalent correlation functions for particles of the same color, 

[(p(x,y)p(0,0))} = hm(p a (x,y)p a (0,0)). (77) 

Can we extract information about this correlation function from the Bose gas 
correlations? Assuming that correlations of the density fluctuations 5p a (x,y) 
have a finite limit for n — > 0, it follows from Eq. (76) that we can expect 
two types of terms in the n — * limit of Eq. (72): Contributions diverging as 
~ 1/n, and those which saturate. Counting powers of n shows that we cannot 
identify the diagonal correlations needed in Eq. (77), but only the difference 
between diagonal and off-diagonal correlations. However, this differene can be 
calculated both in the RG approach and the variational ansatz (VA), and is 
therefore of interest, too. To check for the diverging contributions proportional 
to 1/n in Eq. (72), we have to map the Bose gas exponent rjb to the corre- 
sponding exponent n of the replicated line lattice. Using the result of Eq. (71) 
for the mapped velocity of sound, % — > T and m& — > ng, we get in the limit 
n — > the diverging exponent 




V=— \—c = — yJldP- (78) 



For the last expression, we have used the definition of the length scale Id given 
in Eq. (63). Focusing on the homogeneous part of the line density correla- 
tions, we obtain, by comparing the diverging part proportional to 1/n, the 
correlation function 

[(p a (x,y)(p a (0,0) -M0,0))>] = - P ^4hp ^j^'lyj^y + °*-> ( 79 ) 
where o.t. stands for omitted oscillating terms. The length scales in the last 



expression are a = 1/p, the mean separation of lines, and £y = ygT/p 3 A 
corresponding to 1/ pv s , the mean longitudinal distance between collisions of 
lines. 

Before we discuss the oscillating contributions of the correlation function, we 
would like to compare our result for the difference correlations in Eq. (79) 
with the predictions made by RG techniques and VA. Both, RG and VA are 
constructed on the basis of linear elasticity theory, which describes the line 
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lattice in terms of a scalar displacement field u(x, y). To represent the coupling 
of the random potential to the line density in terms of the displacement field, 
the line density in each replica a has to be decomposed into harmonics as 
[32,65,50] 

oo 

p a (x, y) = p {1 - d x u a (x, y)}J2 cos{2nmp(x - u a (x, y))}. (80) 

m=0 

In this representation, higher order terms in d x u a (x,y) are usually ignored 
since u a (x,y) varies slowly on coarse-grained scales. By retaining only the 
coupling of the first harmonic (m = 1) to the random potential, one obtains 
for the displacement field u(x,y) the Hamiltonian of the RFXY model [42]. 
For this model, the RG approach [44] gives the result 



i(d x u a (x, y)d x u a (0, 0))] = |^ |JL In (r/a) + r 2 In 2 (r/a)} , (81) 

d 2 a 2 



[(d x u a (x, y)d x up(0, 0))] = ^^r 2 In 2 (r/a) , (82) 

where the In 2 contribution is obtained by an expansion to lowest order in 
r = 1 — T/Tg which measures the distance from the glass transition at T g = 
4ttJ . Here J is the elastic constant of the (isotropic) RFXY model and r 2 = 
x 2 + (cn/c^y 2 the rescaled squared distance. J can be expressed in terms of 
the elastic constant Cn and c 44 , which are related to compressions and tilts of 
the line lattice, by [50,55] 

a 2 

J = -^y/C U Cu. (83) 

In contrast, VA [53] yields a simple logarithmic growth of displacements with 
r = 0. The amplitude of the logarithm in Eq. (81) then depends on whether 
one allows for replica symmetry breaking or not. The replica symmetric result 
for the amplitude agrees with the RG result in Eq. (81), as compared to the 
universal amplitude a 2 /2ir 2 , which is obtained in case of replica symmetry 
breaking. 

The relation between density- density correlations and displacement correla- 
tions follows from Eq. (80) as 

[{p a (x, y)p p (0, 0))] = p 2 (1 + [{d x u a (x, y)d xUf} (0, 0))]) + o.t., (84) 

where we have not written explicitly the oscillatory contributions which are 
generated by harmonics of order m > 1 in the density decomposition. Using 
this relation, we can extract displacement correlations from our Bethe ansatz 
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result in Eq. (79). Comparing Eq. (79) and Eq. (84), we deduce 



d 2 a 2 r — / x 2 y 2 ^ 



[(d x u a (x,y)(d x u a (0,0) - 3^(0,0)))] = Q^-^yJldP In \-j + ^ ] (So) 

To complete the comparison with RG and VA, the elastic constants have to be 
determined in the low density limit p <C I /Id, where the above prediction holds. 
The tilt modulus is given by C44 = g/a, whereas the compression modulus Cn 
is determined by the mutual interaction potential W(x) of the lines. Due to 
entropic and steric repulsions, the original contact potential becomes long- 
ranged. The interaction potential can be read off from the free energy density 
f(a = 1/p) since in a grand-canonical description with chemical potential \i 
one has 

f(a) = {-g(n/n c - 1) + W(a)}/a. (86) 

The elastic constant en is now determined by the curvature of the free energy 
density, evaluated at the mean line separation a, as 

c n = x 2 f"{x)\ x=a = aW"(a) = A ( 87 ) 

To calculate the last expression, which is valid in the low density limit a ^> Id, 
we have used the result for the free energy in Eq. (60), which yields the non- 
trivial part W(a)/a of the free energy density. From this result for the elastic 
constants we obtain J = (agA/T) 1 / 2 / An 2 , and hence the amplitude of the In 
term of the RG result in Eq. (81) becomes 

a 2 T a 2 r— 

= (88) 

which is in complete agreement with the Bethe ansatz result of Eq. (85). More- 
over, there is also agreement between the two results regarding the anisotropy 
parameter, since £y = (gT '/ 'p 3 A) 1 / 2 = a(c 44 /cn) 1 / 2 in the dilute limit. This 
means that the velocity of sound is related to the elastic constants by the 
mapping between bosons and replicated lines according to 

cn X 1 / 2 

— <— v s . (89) 
c 44 / 

It is interesting to note that exactly the same relation holds in the case of a 
pure 1 + 1-dimensional line lattice, which can be mapped to world-lines of free 
fermions [14] with v s = nhp/m <-> nTp/g. 

Compared to the VA, our Bethe ansatz result for the difference correlations 
[(d x u a (x, y) (d x u a (0, 0) — d x up(0,0)))] is only in agreement with the replica 
symmetric case. A universal amplitude of the logarithm as predicted by the 
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VA with replica symmetry breaking is not reproduced by our approach. How- 
ever, the agreement of the Bethe ansatz result with the RG predictions is not 
complete. If the RG is assumed to give the correct result, then there should 
appear also a slower decaying contribution proportional to ln(r)/r 2 in the 
Bose gas correlations in Eq. (72) which remains finite in the n — > limit. This 
contribution corresponds to the last off-diagonal term in Eq. (76), for which 
this slower deacy follows from the RG result in Eq. (82). 

There are two potential reasons for the absence of such kind of contributions 
in our Bethe ansatz approach. The first one is not related to our mapping, but 
can be traced back to the fact that there are indeed logarithmic corrections 
to Haldane's harmonic fluid approach for the ID Bose gas. These corrections 
vanish in the strong coupling limit q, — > oo and have been calculated explic- 
itly so far only in a 1 /c& expansion by Korepin [66] using the quantum inverse 
scattering method and by Berkovich and Murthy [64,67] by means of confor- 
mal field theory. However, these results predict, in agreement with each other 
to first order in 1/q,, logarithmic corrections only to the oscillating term of 
Haldane's result in Eq. (72). Therefore these results do not account for In con- 
tributions to the homogeneous part of the density- density correlation function. 
But we are interested in the weak coupling limit c& ~ n 2 — > 0, where the actual 
form of the corrections is not yet known. The second potential, and we think 
more likely reason for the absence of corrections is related to the fact that 
we apply our mapping to obtain the excitation spectra of the U(n) fermions 
from the Bose gas spectra by assuming that only non cluster breaking replica 
symmetric excitations are gap- less. This assumption is certainly justified for 
integer n due to the finite energy necessary to break up a cluster. But if n 
goes to zero, cluster breaking excitations might become gap-less, providing a 
possible source of additional contributions to the asymptotic density-density 
correlations. Note that these cluster breaking excitations can be absolutely 
replica symmetric since n particles of mutual different colors can be trans- 
fered from different clusters across the Fermi surface to form together a color 
neutral excitation. Taking into account this type of excitations is, of course, 
very interesting^], but not possible within the mapping to bosons, which itself 
cannot be broken apart. 

Finally, we come back to the oscillating term in the Bose gas correlation func- 
tion of Eq. (72), and its implications for the correlations in the line lattice. 
This term and higher harmonics are quite important in view of experimental 
determinations of the translational order of line lattices. Scattering experi- 
ments measure the Fourier transform of the density- density correlation func- 
tion around a reciprocal lattice vector 2irmp with integer m. Therefore, these 
experiments can determine the value of the line lattice analogue rj, of the ex- 



9 The splitting of clusters has been considered by Parisi for a single line [N = 1) 
to study the possibility of replica symmetry breaking [68] . 
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ponent rjb of the Bose gas correlations. Using the mapping between lines and 
bosons, we have seen that i] b ~ 1/n [see Eq. (78)] and therefore 

77 -> 00, (90) 

in the n — > limit. This result is in strong contrast to the case of a pure 
line lattice with rj = 2. The infinite exponent means that the oscillating terms 
in the density correlation function decay faster than any power law due to 
destruction of quasi-long-range order by the random potential. Exactly the 
same result was obtained by Cardy and Ostlund in their RG approach [43], 
but it is in contradiction to all VA approaches. Although a mean squared 
relative displacement which grows slightly faster than a logarithmic increase 
cannot be obtained directly from the homogeneous part of the Bose gas density 
correlations, the fact that the Bose gas exponent rjb is mapped onto an infinite 
rj for the line lattice might be interpreted as a signature of the presence of 
these contributions. However, this argument has to be viewed with care since 
77 ~ 1/n means that the oscillating terms of the density-density correlation 
function have an essential singularity in the n — > limit. 



7 Conclusion and discussion 

In this paper, we studied the probability distribution function of the free en- 
ergy of a lattice of lines, distorted due to thermal and disorder fluctuations. Us- 
ing the replica approach in combination with an analytically continued nested 
Bethe ansatz, we calculate exactly the cumulants of the free energy in the 
dilute limit. Our results confirm a previously proposed scaling form for the 
replica free energy in terms of the the replica number n, as a scaling field [17]. 
These exact results were obtained by establishing a novel analogy between 
the random line lattice and the weakly interacting ID Bose gas with delta- 
function interactions. Based on this analogy, we related the results for the 
density-density correlation function of the Bose gas, known from Haldane's 
effective theory of ID quantum liquids, and from conformal field theory, to 
correlation functions of the line lattice. Our result for the difference between 
diagonal and off-diagonal replica-correlations is in agreement with RG predic- 
tions, and a replica symmetric variational ansatz for an elastic model of the 
line lattice, if one takes into account the correct disorder dependence of elas- 
tic constants obtained by the Bethe ansatz approach for the free energy. The 
Bethe ansatz results do not agree with the variational ansatz with replica sym- 
metry breaking. The comparison with the RG results of Cardy and Ostlund 
[43] and Goldschmidt and Houghton [44], which predict a faster than log- 
arithmic growth of relative line displacements, shows that our mapping to a 
Bose gas cannot capture the corresponding contributions. Potential sources for 
these additional logarithmic terms within the Bethe ansatz provide avenues 
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for future explorations. 



Finally, it might be interesting to note that all results obtained in this paper 
are valid if the length scale lj, setting the crossover to random behavior for 
a single line, is much larger than the correlation length £o of the random 
potential. When this condition is fulfilled, £o does not show up in the final 
results, and can be safely set to zero from the beginning as we did throughout 
the paper, see Eq. (3). However, if temperature is decreased, the length scale 
l d = T 3 /gA also decreases, and gets equal to £ at the crossover temperature 
T* = (^ofi'A) 1 / 3 , which is of the order of the smallest random energy barrier 
[69]. Therefore, our results apply directly to the high temperature regime 
T ^> T*. They can also be extended to the low temperature limit T -C T*, 
using a method developed by Korshunov and Dotsenko [70] to calculate the 
replica free energy of a single line in a random potential with finite £ . By 
splitting the n-cluster for a single line into n/k separate blocks with k ~ T/T*, 
they were able to show that the ground state energy of their model with finite 
£o can be obtained from the model with 5-function interactions by the simple 
substitutions g — > kg, A — > A; 2 A and n — > n/k. Since the internal structure 
of individual clusters remains unchanged by considering a finite density of 
clusters (see the ansatze in Eqs. (30)-(32)), our results can be translated easily 
to the low-temperature regime by the above substitutions. 
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